A v = 2/5 Paired Wavefunction 
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We construct a wavefunction, generalizing the well known Moore-Read Pfaffian, that describes 
spinless electrons at filling fraction v — 2/5 (or bosons at filling fraction v = 2/3) as the ground state 
of a very simple three body potential. We find, analogous to the Pfaffian, that when quasiholes are 
added there is a ground state degeneracy which can be identified as zero-modes of the quasiholes. The 
zero-modes are identified as having semionic statistics. We write this wavefunction as a correlator 
of the Virasoro minimal model conformal field theory M(5, 3). Since this model is non-unitary, we 
conclude that this wavefunction is a quantum critical state. Nonetheless, we find that the overlaps 
of this wavefunction with exact diagonalizations in the lowest and first excited Landau level are very 
high, suggesting that this wavefunction may have experimental relevance for some transition that 
may occur in that regime. 



I. INTRODUCTION 

The vast majority of quantum Hall states observed 
experimentally in the lowest Landau level (LLL) are 
very accurately described in terms of composite fermion 1 
(or equivalently the hierarchy 2 ) wavcfunctions. Despite 
these successes, there are a number of quantum Hall 
states that remain much less well understood and may 
require more exotic explanations. For example, there ex- 
ist a few experimentally observed quantum Hall plateaus 
in the LLL that do not fit into the usual hierarchy- 
composite-fermion framework 3 . Also, in the first excited 
Landau level (ILL), even the "simple" observed filling 
fractions (such as v = 2 + 1/3 and 2 + 2/5) appear from 
numerics to have strong differences from the correspond- 
ing states in the Of the non-hierarchy exotic 
states that have been proposed, perhaps the best under- 
stood is the Moore- Read Pfaffiar. 6 , which is thought to 
describe^ the plateau observed at v = 2 + 1/2, and whose 
quasiparticle excitations have exotic nonabelian statis- 
tics. However, for the neighboring experimentally ob- 
served plateau^ at v = 2 + 2/5 there are at least two com- 
peting trial states which have been proposed: the hierar- 
chy (composite fermion) stated, and the (particle-hole 
conjugate of the) Z3 Read-Rezayi parafermion stat o?' 1 ? — 
a generalization of the Pfaffian which has an even richer 
nonabelian structure. Another case where more exotic 
states may occur is in the quantum Hall effect of rotat- 
ing bosonsii. 

In this paper we will study another type of generaliza- 
tion of the Pfaffian that gives a different trial state at 
v = 2/5 (or 2 + 2/5) for spinless electrons. We call this 
new wavefunction the "Gaffnian" , for reasons described 
below. The Hamiltonian that generates this state as its 
unique highest density zero energy state is an extremely 
simple generalization of the Hamiltonian that similarly 
generates the Moore-Read Pfaffian (and is also similar in 
spirit to the Hamiltonians that generate the Read-Rezayi 
wavefunctions) . Similar to the Pfaffian and Read-Rezayi 



states, when additional flux is added, there is a degen- 
eracy of states associated with zero-modes of the quasi- 
holes. In the present case the zero-modes have semionic 
statistics, compared to fermionic statistics in the Pfaffian 
case, or parafermionic statistics in the Read-Rezayi case. 
We then write this wavefunction as the correlator of a 
conformal field theory, known as the Virasoro A4(5, 3) 
minimal model. Since this field theory is nonunitary 
we conclude that the wavefunction does not represent 
a phase, but rather a quantum critical point between 
phases. As such, the concept of nonabelian statistics 
is not necessarily applicable. Nonetheless, exact diag- 
onalizations show extremely high overlaps with this trial 
wavefunction. We take this to suggest that this wave- 
function is likely relevant to a phase transition that is 
somehow "near" the hierarchy phase. 

The outline of this paper is as follows. In sectionlTTlwe 
will introduce this new wavefunction and explore some 
of its properties. We define this state as the unique 
highest density zero energy state of a simple Hamilto- 
nian. We then consider what happens when additional 
flux is added to the system. As mentioned above, in the 
presence of quasiholes, there is a ground state degener- 
acy stemming from semionic zero-modes. Since some of 
the analytic manipulations are messy, we relegate some 
of the details to rather lengthy appendices. In partic- 
ular, the demonstration that the Gaffnian is a unique 
ground state of this Hamiltonian, and the explicit count- 
ing of quasihole states is put in appendix^ However, for 
the interested reader, this appendix shows explicitly the 
mechanism by which the semions occur. In section ITTT1 we 
construct the Gaffnian as a correlator of the M(5, 3) Vi- 
rasoro minimal model conformal field theory. In section 
II VI we examine results of exact diagonalizations. We look 
at low energy excitations to find evidence of criticality, 
and we also discover that the overlap of the Gaffnian with 
the hierarchy wavefunction is remarkably high (we also 
find high overlap with exact diagonalizations of systems 
with interactions close to Coulomb). Finally, in section 
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[V]we give a brief discussion of some of our results. 

II. THE GAFFNIAN WAVEFUNCTION 

Before constructing our new trial wavefunction, for 
motivation we review construction of Laughlin wavefunc- 
tions. Constrained to a single Landau level (LLL or ILL), 
the relative angular momentum of two fermions L>2 must 
be odd and positive. Thus the minimum relative angular 
momentum is L 2 nm = 1 (For bosons, L™ n = and L2 
must be even). We define a projection operator Pf to 
project out (i.e, to keep only) states where any two par- 
ticles have relative angular momentum less than L 2 nm +p 
(with p even). This projection operator— can serve as a 
Hamiltonian. The Laughlin v = \j(L 2 nm +p) state is the 
unique highest density (zero energy) ground state of Pf 
(with L™ m + p odd for fermions and even for bosons) . 

We now generalize this construction. In a single 
Landau level the relative angular momentum of three 
fermions L3 has a minimum value L™ in = 3 (For 
bosons, L™ in = 0). It is not hard to show 14 ! 15 . 16 that 
L 3 7^ L™ in + 1 is dictated by symmetry, but any other 
L3 > L™ m is allowed. We analogously define a projec- 
tion operator P3 to project out (i.e., to keep only) states 
where any three particles have relative angular momen- 
tum less than L™ m + p which will serve as our Hamilto- 
nian. It is well known 6 i 13 i 14 that the Pfaffian (at v=l/2 
for fermions and v = 1 for bosons) is the unique highest 
density (zero energy) ground state of the Hamiltonian 
P3 2 . In Ref. ITEI another state, known as the "Haffnian" 
is shown to be the unique highest density (zero energy) 
ground state of P3 (which is a v = 1/3 state for fermions 
and v = 1/2 state for bosons). Using the method of 
we can show that the Hamiltonian Pf also has a 
unique highest density (zero energy) ground state. The 
argument is straightforward and is given in the appendix 
of this paper. This unique state occurs at v = 2/5 for 
fermions (y = 2/3 for bosons), and is the focus of this 
paper. Since this new p = 3 state lies between the p = 2 
PFaffian and the p = 4 Haffnian, we alpha-phonetically 
interpolate and dub this new state the "Gaffnian" . 

Before commencing our study of the Gaffnian, we note 



that several other states can be constructed analogously. 
By considering general fc-particle angular momenta Lk, 
we can construct a general Pf . In Ref. it was shown 
that the Hamiltonian P| generates the Zk-i Read-Rezayi 
state (the Z2 state being the Pfaffian). Study of several 
other values of p and fc is given in Ref. [16| by three of the 
current authors. 

Knowing that a unique quantum Hall ground state ex- 
ists for the Hamiltonian P3 , we set about describing the 
properties of the Gaffnian. We begin by writing down 
the wavefunction explicitly. 

We will represent a particle's coordinate as an ana- 
lytic variable z — x + iy which is simply the complex 
representation of the particle position r. On the spher- 
ical geometry, z is the stereographic projection of the 
position on the sphere of radius R to the plane. All dis- 
tances are measured in units of the magnetic length. We 
can write any single particle wavefunction as an analytic 
function times a measure jtx(r). On the disk the measure 
is the usual gaussian factor- /i(r) = e~l z / 4 whereas on 
the spheroid (with stereographic projection to the plane) 
the measure is fi(r) = [1 + |z| 2 /(4P 2 )r (1+JV< * /2) with 7V 
(=2P 2 when the magnetic length is unity) being the total 
number of flux penetrating the spher o 9 ! 16 . On the sphere 
the degree of the polynomial ip(z) ranges from z° to 2;^ 
giving a complete basis of the + 1 states of the LLL. 
We will not write the measure explicitly, instead writ- 
ing all wavefunctions simply as analytic functions (which 
must be fully symmetric for bosons and fully antisym- 
metric for fermions). 

It is convenient to think, for a moment about bosons at 
v = 2/3. Since the Hamiltonian P3 puts no restriction on 
the two particle angular momentum, there is no restric- 
tion against two bosons being at the same point zq. How- 
ever, when a third particle approaches, it must approach 
the other twoi£ such that the overall angular momen- 
tum of the three particles is p > 3, i.e., the wavefunction 
vanishes as (23 — z$) p . (In this sense, the Gaffnian, like 
the Pfaffian and Haffnian, is a paired state in the spirit 
of that originally proposed in Ref. 0.) The Gaffnian 
wavefunction can be written explicitly as 18 
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where q — corresponds to a bosonic (v — 2/3) wave- 
function and q = 1 is a fermionic {y = 2/5) wavefunction. 
We have assumed the number N of particles is even and 
S means symmetrize or antisymmetrize over all particle 
coordinates for bosons or fermions respectively. One can 



confirm directly that the above wavefunction for q = 
correctly has the property that it does not vanish as two 
particles come to the same position but vanishes as three 
powers as the third particle arrives. Further, we show 
in appendix that this is the unique densest wavefunc- 
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tion that has this property. As is standard in spherical 
geometry 2 *^ we can obtain the value of the flux by look- 
ing at the maximum power of Zi that occurs. Counting 
powers of z we find that the Gaffnian wavefunction oc- 
curs on the sphere for total flux 



N^=3N/2-3 + q(N -1). 



(2) 



This value of flux is the same as that for the standard 
hierarchy v = 2/5 stated. This should not be a surprise, 
since some of the first trial wavefunctions for v = 2/5 
(for fermions) in the hierarchy were based on pairing 17 . 
In the appendix we analytically establish that this state 
is the unique zero energy state of the Hamiltonian Pg 
at this flux. We have numerically confirmed this fact by 
explicitly diagonalizing the Hamiltonian Pg with up to 
N = 12 particles on the spherical geometry and up to 
N = 10 particles on the torus. (We note that on the 
torus the Gaffnian occurs at = (3/2 + q)N meaning 
there is no "shift" , which is always the case on the torus) . 

In appendix iGl we consider possible generalizations of 
the form of the Gaffnian wavefunction (Eq. ^) . In par- 
ticular, we find trial states for wavefunctions of the Jain 
series v = p/ (mp+1) (with m odd for bosons and even for 
fermions) with the same value of the flux as the usual Jain 
sequence. Since (as we will discuss below) the Gaffnian 
is distinct from the hierarchy (or Jain) states, we sus- 
pect that these trial states are similarly distinct from the 
usual Jain states. However, we leave detailed study of 
these wavefunctions for further work. 

Since the Gaffnian is a paired Btat oPi , we expect that 
each additional flux added will correspond to two quasi- 
holes, each with charge e* = ev/2 with — e the charge 
on the elementary underlying "electron" (or underlying 
boson for q = 0). Generally, we define the number of 
extra flux added to the Gaffnian ground state to be 



N$ — (3N/2 — 3 + q(N — 1)) 



(3) 



(compare Eq. [5J. Note that n here is defined so that it 
is half integer if N is odd. To construct wavefunctions in 
the presence of n (integral) additional flux, we can insert 
a factor of 

n ( z *- w j) n Ofc-w fe ) (4) 

a<N/2;j<n N/2< b<N :n<k<2n 

into the above wavefunction inside the symmetrization 
where the u>'s indicate the quasihole positions. However, 
for 2n fixed quasihole positions, there are apparently ( 2 ^ 1 ) 
inequivalent ways to choose which of the positions Wj are 
labelled with an index j < n and which with an index 
j > n. One might expect that the different groupings of 
the positions into these two groups generate equally many 
inequivalent quasihole wavefunctions. The fact that we 
find more than one independent quasihole wavefunction 
means that there are zero-modes associated with these 
quasiholes^. Analogous to the Pfafnar. 14 >, however, it 
turns out that there are many linear dependencies be- 
tween these many different wavefunctions. 



In appendix ^ we show explicitly how to count the 
zero energy ground state degeneracy of the system with 
Hamiltonian Pg at any flux (Strictly speaking the ap- 
pendix only addresses the case of N and 2n even, al- 
though the odd case proceeds similarly). We find that 
the degeneracy of zero energy states is given by 



F max 

E 

-i) f =(-i) 



F)/2- 
2n 



2n 



F/2-1 
F 



(5) 



where the maximum value of F is given by F max = 
min(iV, 2n— 2). To verify this result, we have numerically 
performed exact diagonalizations. For every case we have 
examined, we find perfect agreement between this ana- 
lytic rule and the results of our exact diagonalization of 
the Hamiltonian Pg. (We have examined N — 4, 6 with 
n < 6, N = 8 with n < 4, N = 10 with n < 2, N = 5 
with n < 7/2 and N = 7 with n < 5/2.) 

The first term in Eq. [S] corresponds to the positional 
degeneracy of the quasiholes and can be thought of as 2n 
bosons in (N — F)/2 + 1 orbitals. The second term is 
the degeneracy of the zero-modes and can be thought of 
as F fermions in n + F/2 — 1 orbitals. Since the number 
of orbitals changes half as fast as number of particles, 
these zero-modes are a realization of semionic exclusion 
statistics^. 

The form of Eq.[S]is quite analogous to the zero-mode 
counting expressions found for the Pfaffian 14 , Haffniani^, 
and Read-Rezayi State s 9 ! 20 . However, in those cases the 
zero-modes have fermionic, bosonic, and parafermionic 
statistics respectively. For the fermionic (Pfaffian) case 
we put F fermions in a fixed number n orbitals^. For 
the bosonic (Haffnian) case^, we put F fermions in n + 
F — 2 orbitals (which is equivalent to putting F bosons 
in n — 1 orbitals). The Gaffnian case is quite naturally an 
interpolation between these two cases. (The Read-Rezayi 
parafermion case cannot be phrased in this language so 
easily). 

As with the Pfaffian 14 , Haffnian 15 , and Read- 
Rezayi 9 ^ cases, the structure of Eq. |S]also tells us how 
to decompose these degenerate states into angular mo- 
mentum multiplets. We simply calculate the multiplets 
of the 2n bosons in (N — F)/2 + 1 orbitals and also the 
multiplets of the F fermions in n + F/2 — 1 orbitals and 
then add these together using standard angular momen- 
tum addition rules. An explicit example of this angular 
momentum addition is given in appendix iBl 

As discussed above, we can also look at wavefunctions 
with fixed quasiparticle positions. The number of linearly 
independent states should just be given by the zero-mode 
contribution to the above equation 



F,(-i) F =(~iy 



n + F/2 - 1 
F 



(6) 



Indeed by generating wavefunctions (described by Eq. 0] 
inserted into Eq. ^) numerically and checking for linear 
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independence, we find precisely this number of indepen- 
dent states for all cases we have tried (N = 4, n < 6; 
N = 6,n < 5, N = 8,n < 3). 

It is interesting to note that in the case of N > 2n — 2 
(so F max — 2n — 2) the sum Eq. gives the 2n — 1 st 
Fibonacci number, Fib(2n - 1) = Fib(N qh - 1). This 
can be proven trivially by induction on n to show that 
D n + D n+ i/2 = D n+ i. We note that the Z$ Read-Rezayi 
state where also has a degeneracy of Fib(N q h — 1)- An- 
other similarity we have found is that both states have a 
2-fold degeneracy of the ground state at zero momentum 
on the torus geometry (in addition to the usual center 
of mass degeneracy^). However, the two ground states 
occur at different values of the flux for a finite spherical 
system, so they are topologically different states. Also, 
as mentioned above, the state counting formula anal- 
gous to Eq. [S] involves parafermionicSS zero-modes for 
the Read-Rezayi case compared to semionic modes for 
the Gaffnian. 





A 




3/4 


f 


1/5 


a 


-1/20 



X 


<P 


V 


a 





1 








a 


1 + if 




a 


<P 


a + V> 


1 + <p 



FIG. 1: In the Virasoro minimal model conformal field theory 
-M(3, 5), there are three nontrivial fields, ip, ip and a with 
dimensions A given in the left table and fusion algebra given 
in right table. 



III. CONFORMAL FIELD THEORY 

We now write this Gaffnian wavefunction as a correla- 
tor of a conformal field theory (CFT)Si. Making the con- 
nection to CFT has, in the past, been extremely powerful 
in understanding states with nonabelian statistics (See 
for example Refs. IHHH. For example, the structure 
of a CFT can tell us about behavior of the degenerate 
space under adiabatic braiding of quasiholes 2 ^. A CFT 
describing a paired state should contain a field ip with fu- 
sion relation ip x ip *~ 1 such that it has operator product 
expansion 



ip(z)ip(w) ~ (z — w) 



-2A„ 



[1 



(7) 



with 1 the identity, A^, the conformal weight (or dimen- 
sion) of ip, and dots representing less singular terms. We 
can then construct a paired wavefunction 



*=(n^))n^ 



z,) 2A * +9 



(8) 



Repeating the arguments which are presented in Ref. it 
is clear that (for q = 0, i.e., for bosons) this wavefunction 
will not vanish as 2 particles come to the same position 



since the (fractional) Jastrow factor precisely cancels the 
singularity of the operator product expansion. However, 
the wavefunction vanishes as z 4A * powers when the third 
particle approaches the other two (since there are three 
(fractional) Jastrow factors and only one singularity). 
The Moore- Read PfaffianS is described in this way by the 
Ising CFT, also known as the A4(4,3) minimal modelSi, 
which contains such a field ip with weight A^ = 1/2 so 
the wavefunction vanishes as z 2 as three particles come 
to the same point. The Gaffnian is correspondingly de- 
scribed by one of the simplest generalizations of the Ising 
CFT, known as the minimal model At (5, 3). This CFT 
has a field ip with A^ = 3/4 so that the wavefunction 
vanishes as z 3 as three particles coalesce (for q = 0). 
The dimensions and fusion rules for the three indepen- 
dent fields (ip, a and ip) in this model are given in Fig. ^ 
The fusion of the field a with the field ip gives us the 
operator product expansion^! 



ip(z)a(w) ~ (z — w) 1 / 2 ip(w) + 



(9) 



where here the exponent —1/2 is determined by the con- 
formal weights in Fig. [3 as A v — A^, — A CT . As described 
in Ref. this power of 1 /2 means that the quasihole cre- 
ated by the field a must have charge e* = ev/2 consistent 
with our expectation for a paired state. To see how this 
happens we write a general wavefunction in the presence 
of 2n quasiholes as 



N 



i<3 



i=l I 

N 2n 

> 3/2+9 nn^-^) i/2 (io) 

i=li=l 



Given the operator product expansion Eq. [51 the final 
exponent in Eq. [!|] must have power 1/2 so that the 
wavefunction is single valued in the z's. This Jastrow 
factor then pushes precisely a charge ev/2 away from 
each quasihole. 

We can also use the fusion rules to count the degener- 
acy of the 2n quasihole state. The degeneracy is given 
by the number of ways the a fields in Eq. EI] can fuse to- 
gether to form the identity. This is illustrated graphically 
as the number of paths through the Bratteli diagram 22 
shown in Fig. [21 The number of paths is Fib(2n — 1), 
which is consistent with the result of our above counting 
formula. If the number of particles iV is even, then we 
pair the ip fields to form identities, and the a fields must 
also pair to form the identity. However, if N is odd, we 
can only form the identity if the a fields fuse to form one 
more ip that can then fuse to form the identity with the 
one remaining ip field. 

One may ask how we know that we have the correct 
conformal field theory (particularly in light of the fact 
that classification of all conformal field theories is an on- 
going research field). The fact that we have a paired 
state at filling fraction v — 2/3 for bosons (i.e., the fact 
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FIG. 2: The Bratteli diagram shows how the 2n quasihole 
fields a fuse together. This is just a graphical representation 
of the fusion rules (Table ^ where at each horizontal step, 
the states at the previous horizontal position are fused with 
one more a field. The number of conformal blocks — which 
gives the nonabelian degeneracy — is seen graphically by the 
number of paths through the diagram starting and ending at 
the bottom when N (and 2n) is even. When TV (and 2n) is 
odd, the path needs to start at the bottom but end at the 
top to fuse with the one unpaired t[> field. By straightforward 
counting, the number of such paths with 2n steps can be seen 
to be the 2n — 1 st Fibonacci number. 



that the wavefunction does not vanish when two particles 
come together) means we must have a field ip which fuses 
with itself to form the identity. The fact that the Hamil- 
tonian forces the wavefunction to vanish as three pow- 
ers when three particles come together further fixes the 
dimension A^. It is easy enough to show that the only 
Virasoro minimal model conformal field theory with such 
a field is M(5, 3). If we further insist that the charge of 
the quasihole should be ev/2, as is expected for a paired 
state, this fixes the exponent of the final factor in Eq. 
1101 and this in turn fixes — A CT = 1/4. We must also 
insist that the fusion relations for fusing many quasipar- 
ticles with each other has the form of the Bratteli dia- 
gram in Fig. [21 Finally, one can look at the subleading 
behavior of the wavefunction as particles approach each 
other to extract the central charge of the theory, which 
again is consistent with Al(5, 3) (we do not perform this 
calculation here). These restrictions place serious con- 
straint on any possible conformal field theory we would 
like to use to represent the gaffnian state. Certainly there 
is no "simple" (i.e., minimal model) theory other than 
M(5, 3) with the required properties. However, we have 
not proven that no other theory exists. 

The conformal field theory A4(5,3) is nonunitarySi. 
This highly suggests that the Gaffnian wavefunction does 
not represent a true phase, but rather represents a quan- 
tum critical point. The argument for this goes as follows: 
The edge state theory in 1+1 dimensions of a quantum 
Hall state should be described by the same conformal 
field theory as the bulk 2 dimensional theory. However, 
since the edge state theory is a dynamical theory, it must 
be unitary. If we have a trial wavefunction that is gen- 
erated by a nonunitary theory, apparently the only way 
out of this conundrum is that the edge state theory does 
not exist; i.e, edge excitations do not stay on the edge, 
but leak into the bulk. This could indeed be the case 
if the ground state has arbitrarily low energy excita- 
tions in the thermodynamic limit. This could in turn 



occur if the wavefunction represents a quantum critical 
point. Indeed, there have been past examples of critical 
quantum Hall states which are described by nonunitary 
CFTs 14 i 15 i 24 . While there is no strict proof that a nonuni- 
tary conformal field theory necessarily implies a critical 
state, there is also no understanding of how anything else 
could occur. 



IV. EXACT DIAGONALIZATIONS 

We now turn to exact diagonalizations. Strictly speak- 
ing, the Hamiltonian P| has been defined to be a projec- 
tion operator that acts on the full wavefunction (to keep 
any states where any three particles have relative angular 
momentum less than three). As such, this Hamiltonian 
has eigenvalues that are either zero (for the zero energy 
space) or unity. A more physical version of this Hamil- 
tonian can be written as 

H = E ( V 3fiP? jk + V S , 2 P? jk ) (11) 
i<j<k 

where we have defined a general three body operator Pf- k 
which projects out (i.e., keeps) any component of the 
wavefunction where the three particles i, j and k have rel- 
ative angular momentum L m i n + p. (On the spher o 12 ! 16 , 
one defines P[j k to project out (i.e., keep) any cluster of 
three particles with total angular momentum 3N ( / > /2—p). 




i > i i i i i i i i i i i > i 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 

L/N 

FIG. 3: Lowest neutral gap excitation of the Gaffnian as a 
function of system size (using the Hamiltonian in Eq. 1111 in 
units of V3,o = V3,2- Data is shown for N = 8, 10, 12 particles. 
The solid is a linear fit of all three data points. The dashed 
line is a fit of the two larger systems only (suggesting that if 
we could access even larger systems, the extrapolations might 
be even closer to zero). This data suggests the possibility that 
the gap may extrapolate to zero in the thermodynamic limit, 
as would be expected for a critical state. However, from the 
available numerical data, we cannot exclude the possibility 
that it extrapolates to a finite value. 

Note that three particles cannoU^ have relative angu- 
lar momentum of L m i n + 1, so this Hamiltonian gives 
energy to any case where the relative angular momen- 
tum of any cluster of three particles is less than three. ( 
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Some readers may have assumed that the form of Eq. ITT1 
is what we meant all along when we have been writing 
P3 , as we were not very explicit about what we meant.) 
Since the Hamiltonian Eq. II II gives energy to any cluster 
of three particles with relative angular momentum less 
than 3, it has precisely the same zero energy space as 
P3 . However, the excitation spectrum here is different, 
and is dependent on the values of V^o and 1/3,2 • 

Let us first examine the issue of criticality. In Fig. |3 
we show the lowest energy neutral excitation of H (from 
Eq. 111(1 as a function of system size for iV = 8,10,12on 
a spherical geometry with V^o = ^3,2 (We have chosen to 
look at bosons on a sphere because we can go to larger 
systems). As can be seen in the figure it appears that 
the gap extrapolates to a positive value, but it is not 
possible to rule out extrapolation to a zero value which 
would be a sign of criticality. Furthermore, changing the 
ratio of V^o/V^ (data not shown) does not appear to 
substantially affect the ratio of the extrapolated energy 
to the reference energy of the gap for N = 10. 
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FIG. 4: Squared overlaps of trial states with the exact ground 
state at v = 2/5 on a sphere with 10 electrons, as we vary the 
coulomb interaction. Solid line is the overlap of the Gaffhian 
wavefunction with the exact ground state. Lhe dashed line is 
the hierarchy 2/5 state— with the exact ground state. The 
top is results for the lowest Landau level, the bottom is the 
first excited Landau level. In the horizontal direction the in- 
teraction is varied around the Coulomb interaction by adding 
an additional SVi Haldane pseudopotential. 

We now turn to the question of whether the Gaffnian 
is physically relevant to the physics of 2D electron sys- 
tems. We have performed exact diagonalization on a 
spherical geometry for 10 electrons in the lowest Landau 
level (LLL) and first excited Landau Level (ILL), and 
we have varied the electron-electron interaction in the 
neighborhood of the Coulomb interaction by varying the 
Haldane pseudopotential2*i2i coefficient V\. In Fig. 2a, 
we show the overlap of the exact ground state with our 
trial wavefunction. Results are shown for the Gaffnian 
(solid) and the hierarchy 2/5 stat o 2 i 23 (dashed). Over 



a range of Vi both trial states have quite good overlaps 
with the ground state considering that the zero angular 
momentum Hilbert space has 52 dimensions (Note that 
for many of the well known numerical casesi 2 ^ where ex- 
tremely large overlaps have been reported, the dimension 
of the available Hilbert space is much smaller than this). 
Near the regime of V% where the overlaps drop, we be- 
lieve the system is in the Read-Rezayi phased (although 
at a different value of flux on the sphere). Since both 
the Gaffnian and hierarchy states have such large over- 
laps with the ground state, they necessarily have large 
overlaps with each other, although in the thermodynamic 
limit they become orthogonal. 

We have also performed exact diagonalization on the 
torus geometry. Here, the Gaffnian ground state is found 
to be doubly degenerate (in addition to the usual cen- 
ter of mass degeneracy). The two zero energy ground 
states are distinguished by a parity quantum number. 
The state with positive parity again has extremely high 
overlap with the hierarchy state, similar to the overlaps 
on the sphere. As on the sphere, both of these have a 
high overlap with the exact ground state for a wide range 
of interactions. However, we do not find that the exact 
ground state has an even approximate double degener- 
acy in the regimes where the overlaps of the Gaffnian 
and the hierarchy are large. Approximate double degen- 
eracy of the ground state is found where we believe the 
Read-Rezayi state is the proper ground stated. 



V. DISCUSSION 

If the Gaffnian does turn out to be a critical state, as 
suggested here, this then raises the question as to what 
the neighboring phases are. It is reasonable that one 
would be a "strong pairing" phase (albeit one that can- 
not be easily described within BCS theorySi) which may 
correspond to the hierarchy wavefunction itsel£i£. This 
would be quite natural considering the high mutual over- 
laps of the hierarchy and Gaffnian. 

The nature of the state on the opposite side of the tran- 
sition is a bit harder to guess at. One possibility is that 
it is the Read-Rezayi state. This would make some sense 
because of the similar ground state degeneracy. Here, we 
imagine that as we approach the transition from the hi- 
erarchy side, the putative zero energy state would drop 
continuously and hit zero energy at the Gaffnian criti- 
cal point. It would then stay at zero energy through the 
Read-Rezayi phase. On the other hand, we should note 
that there is a notable topological difference between the 
Read-Rezayi and Gaffnian state, which is more evident 
on the sphere as they occur at different values of flux. 

Yet another possible candidate for a state that might 
occur nearby is a charge density wave state. We leave 
the project of sorting out the details of this transition for 
future work. 

If the Gaffnian is in fact a critical state, this means 
that the concept of "nonabelian" statistics^ may not be 
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well defined. Indeed, the idea of statistics describes what 
happens to a system when particles are adiabatically ex- 
changed. Since definition of adiabatic usually requires 
any perturbation to the system to be on a timescale 
slower than h/A with A the minimum gap in the system, 
if the system has gapless excitations, there is generally 
no way to have adiabaticity. One might ask whether any 
remnant of the idea of nonabelian statistics still remain. 
This is a question that is hard to answer without knowing 
the details of what these "critical" low energy excitations 
are. 
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APPENDIX A: ANALYTIC COUNTING OF 
ZERO ENERGY STATES 

In this appendix we will enumerate all possible zero 
energy states of the Hamiltonian Pg on the sphere with 
any number of particles and at any given flux. Our ap- 
proach will be in two parts. In section [A II we will write 
down a linearly independent set of zero energy states (and 
we will count them). It is this section that shows most 
clearly how the semionic zero-modes arise. Then in sec- 



tion lA 2l we will show that these wavefunctions are indeed 
zero energy states of P| . Finally, we will show in IA 31 
that these states indeed form a complete set of the zero 
energy states. The arguments here are quite similar to 
those given in Refs. ITlm However, here the situation 
is more complicated as our zero-modes are not simple 
fermions and bosons as in those two references. Note 
that throughout this appendix we will focus on the case 
where N (and therefore 2n) is even. The N odd case is 
a relatively simple generalization. 



1. Counting States 

We start with the requirement that ip vanishes as three 
powers as any three particles approach each other. The 
wavefunction Eq. ^ clearly provides one such solution 
(at a given value of flux). We will call this the Gaffnian 
ground state. In this section we propose a more general 
form for wavefunctions when there are some arbitrary 
number 2n quasiholes (or n additional flux) added to the 
ground state, and we will count the number of such states 
that are linearly independent. 

Inspired by the work of Refs. and analogous 

to the Pfaffian, Haffnian and Haldane-Rezayi states, we 
write our proposed wavefunction in a form with broken 
and unbroken pairs. Let us declare that F of the N 
total particles are unpaired. Restrictions on F will be 
determined later. We then propose the following form 
for our wavefunctions: 



n c* 



,2+9 



l<a<fc<- 



n (z c -z d ) 2+ « n (*«-*/) 

% + l<c<d<N l<e<-f</<JV 



1+9 



n 



F N + F 



, W2r 



l<9<- 



-+!' • 



, Z n+f ) 



(Al) 



,1=1 i=JV 



+1 



where as above S either symmetrizes (for bosons, even 
q) or antisymmetrizes (for fermions, odd q) over all of 
the z coordinates. Here we have defined $ to be the 
Read-Rezayi quasihole insertion^: 

$(zi, Z 2 ; Wl, . . . , W2n) = (A2) 

1 " 



mined later. (The sum over r S S2n is the sum over 
permutations of the 2n variables w). We now specialize 
to the case of q = (bosons) for simplicity. Since the 
particles must be at the flux = 3(4^ — l) + nwe will 
deduce that the highest degree of the unpaired particle 
coordinates is n— 1. To see how this is deduced we start 
by defining A = {z%, . . . , zn_ }, and J5 = {z» +1 , . . . , zjq\). 
We then simply count up powers of z/~ appearing in tpc 



and fi is a wavefunction for the zero-modes to be deter- 
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V'g = s[H(z Ai z Aj f H(z Bi z Bj f H(z Ai Z B] )\{ 

i<j i<j i.j i 



<S>(z A% ,z Bi ) 
(z Az -z Bi )l 



l[(z Ai -z Bj ) n(...) 



(A3) 



2(f-l)forA 2(f-l)forB 



for A&B 



both paired & unpaired 



-1+n for A&B 

V 

paired only 



£ for A&B ??? for A&B 



unpaired only 



Here S symmetrizes over all coordinates z. In this equa- 
tion, we have written beneath each term the number of 
powers of zu occurring. Thus adding up the powers, we 
conclude that ft is some polynomial in unpaired coordi- 
nates of degree to^ : < to^ < n — 1 — -j for each unpaired 
coordinate Zi. Notice also, that this puts a restriction on 
F: F < 2n — 2, and since obviously F < N, we obtain 

F < min(2n - 2, N) (A4) 

as written above in the main text. The maximum degree 
of SI occurs when F = and is given by n — 1. 

To see how many linearly independent wavefunctions 
we have for given {N, n, F} we proceed as follows. Wc 
choose N — F (necessarily even here) of the N coordi- 



nates and group them together in pairs {(2^,2:^)} for 
i = 1, ... (TV — F)/2 with a u bi G [1 . . . N] and a t ^ a,-, 
hi ^ bj for i 7^ j and a-i ^ bj for all We then bring 
together the position of the paired paired particles to co- 
ordinates Zi. In other words we set z ai = z\ >i = Zi for i = 
1, . . . (N — F)/2. Taking this limit selects out a particular 
group of terms from the original full symmctrization that 
do not vanish. In particular the nonvanishing terms are 
the terms in which a factor of &(z ai , Zb t ; ■ ■ -)(z ai — Zb^ 1 
appeared for each pair (za^z^). The other terms will 
have a factor of (z ai — z^) in the numerator, and will 
vanish in these limits. After taking these limits we are 
left with 



i> = s' 



Y[{Zi - z 3 f Y[ $(5 fe ; tui, . . . , w 2n ) \\ ( Za ~ z& ) 2 II ( Zc ~ Zd f 



(A5) 



l<a<6<-' 



i + l<c<d<~±£ 



N-F F_ 

2 2 



N-F N+F F_ N+F 

2 2 2 2 

n ri( ze n n (^™~ z /) 3 n n (zg-z h ) 2 n(z 1 ,...,zF-,zN +1 ,...,zN+F) 

1 = 1 e=l m=l /=f+l g=l h= N +1 



where S' symmetrizes over {Si} and 

{z\ , . . . , zf , z n+i , . . . , z n+f } separately. In other 

22 2 

words, 5" is what remains of the original full sym- 
metrization over the N particles. The underlined factors 
contain the dependence of tp on z, and are symmetric in 
{Si}, while the doubly underlined factor is symmetric 



in {z\ zf , z jv+i z n+f \ as well. Thus, the 

22 2 

symmctrization S' reduces to S" which symmetrizes 
over unpaired particles only (because the expression 
is already symmetric in zi), and we can rewrite the 
wavefunction as 



i> = Wlj? S" 



N-F F_ 

2 2 



N-F N+F 

no-.) n n^-^) 3 n ri i*™-*t? 

1 = 1 e=l m=l f—N_ + i 



(A6) 



r 



where V'lj is a Laughlin-Jastrow factor in the unpaired ways be taken to be fully symmetric in its arguments (any 
particle coordinates. We thus discover that fi(. . .) can al- nonsymmetric parts vanish when symmetrized). We can 
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thus think of this as a bosonic wavefunction for the zero- 
modes. However, we've already determined the maximal 
degree of fi(. . .) to be n — 1 — The minimal degree is 
obviously 0, so we have a total of n — -j orbitals in which 

to put F bosons. There are ( ( ""^ +f_1 ) = ( n+ f _1 ) 
such linear independent wavefunctions. This is equiva- 
lent to placing F fermions in n + F/2 — 1 orbitals. Since 
the number of orbitals changes half as fast as the number 
of particles we put in them, these particles have semionic 
exclusion statistics 19 . 



Using L = i(a + (3 + 7) and with N<p = 3(N/2 - 1) we 
obtain in this case L = — 6. 

• Case 2a: i,j £ A, i < j; k £ B with pairing of the 
form (ia)(jb)(ck), i.e. terms of the form 

^{z u Z a ) $(Zj,Z b ) gKfcvffc) 

(z t - Z a ) (Zj - z b ) (z c - z k ) 

Here we have 



2. Zero Energy 

We will continue on to demonstrate that the linearly 
independent set of wavefunctions we have just written 
down is in fact a complete set of zero energy states of 
the Hamiltonian P3 . First, however, we show that these 
wavefunctions are indeed zero energy states. The wave- 
function for any zero energy state must vanish as three 
or more powers when three particle positions come to the 
same point. On the spher o 12 ! 16 , this is equivalent to re- 
stricting the total angular momentum of the cluster of 
three particles to be no greater than ZN$/2 — 3. 

First we'll show that the proposed wavefunctions ipc 
are zero energy eigenstates of this Hamiltonian. For the 
ground state, i.e. no additional flux [n = 0) we have 
A^ = 3(-j — !)• Let us look at the (ijk) = [zi,Zj,Zk) 
triplet. We want to know what the highest total angu- 
lar momentum is for this triplet in our wavefunction ipQ. 
The wavefunction can be rewritten (a la HaldaneiS) as a 
sum of terms proportional to / re i(^, Zj, z k )f tot (zi, Zj, z k ) 
where / ro i() is an eigenstate of Uj k , the 3 particle relative 
angular momentum operator, and /tot() is an eigenstate 
of Lijk, the 3 particle total angular momentum operator. 
Note that, as above, we will always focus on q = for 
simplicity (The q =/= case is a relatively minor general- 
ization). To find the total angular momentum, we look at 



a = deg 2i = 2[(- - 1) - 1] + (— - 1) - 1 
(3 = deg z . tJjq = 2[y - 2] + - 1) - 1 (A9) 



N 



7 = deg Zk = 2[- - 1] + (- - 2) - 1 



Similarly, adding up these powers we obtain an angular 
momentum L = |A^j — 4. 

• Case 2b: i,j £ A, i < j; k e B with pairing of the 
form (ia)(jk), i.e, terms of the form 



$(zi,Za) $(z 3 ,z k ) 
(Zi - z a ) (Zj - Zk) 



(A10) 



Here we have 



N N 
a = deg Zi ipQ = 2[(- - 1) - 1] + (— - 1) - 1 

P = deg^Vc = 2[y - 2] + (y - 1) (All) 



N 



7 = deg zt Vg = 2[- - 1] + (— - 2) 



which results in an angular momentum L = jN^ — 3. 



the maximal degree of z^zjzj in /tot (), and find the total These cases are the only possibilities. Thus the highest 



angular momentum L — max(i(a + (3 + 7)). To find the 
maximum total angular momentum we must consider all 
possible ways to have chosen the triplet (zi,Zj, Zk) from 
the many terms in the wavefunction. In particular, we 
must look at all cases of which coordinate is one of the 
paired variables, and which is unpaired, as well as look- 
ing at which variable is an ^4-coordinate, and which is 
a B coordinate. Here we are looking at the relative an- 
gular momentum of a given triplet in each of the many 
terms of the symmetrization sum. All possibilities are 
enumerated next. 



• Case 1: k £ A, i < j < k. 






Here we have 






N 

a = deg z . ip G = 2[(y - 1) 


-2]- 




(3 = deg z . Vg - 2[(y - 2) 


- 1] - 




N 

7 = deg Zfc Vg = 2[y - 3] - 


N 
' ~2 ' 


- 1 



total angular momentum for any triplet is ^N^ — 3 and so 
the proposed wavefunction ipc is a zero energy eigenstate 
of the Gaffnian Hamiltonian P3 as claimed. 



3. Completeness 

Now wc show that the proposed wavefunctions span 
the complete set of zero energy states of the Gaffnian 
Hamiltonian. To do this we will construct the most gen- 
eral zero energy eigenstate and show that it takes the 
form of our proposed wavefunction. Take the follow- 
ing zero energy wavefunction ipc — Wlj^g where here 
ipLj is the Laughlin-Jastrow factor for all of the parti- 
cles. Consider the behavior of tpG as particles in an arbi- 
trary triplet (ijk) approach each other, while the other 
particles remain far away from the three. 
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tp G oc (z t - Zjfizj - z k f{z k - z,f (zi - Zj) q ^ { Zj - z k ) q ^ {z k - Zi) qki (A12) 



a+/3+7=6, part of i))' 2 ., a+l3+l=qij+qjk+qki=Q, part of <j>G 



The wavefunction vanishes as 6 + Q powers as these three 
particles come together. This is equivalenti2ii& to saying 
that the total angular momentum of three particles is L = 
|iV0 — (6 + Q) (Since we're on the sphere, the maximum 
angular momentum of each particle is Any relative 
angular momentum reduces the total by a corresponding 
amoun1pA2ii&) . Furthermore, by analyticity ofipc we must 
have q mn > —2. 

Now, in order for ipG to be a zero energy state of the 
Gaffnian Hamiltonian, we must have Q > —3 (so that the 
relative angular momentum of the cluster is greater than 
or equal to 3 = 6 + Q). From here on we'll concentrate on 
the 4>c factor of the eigenstates, restoring the ubiquitous 
at the end. Allowed forms in the Laurent expansion 
of <Pg as (ijk) approach each other are 



1 



(Zi - Zj) 2 (zj - z k ) 

Zk ~ Zj 

(Zi - z J ) 2 (z j - z k ) 2 
1 

(Zi - Zj)(zj - z k )(z k - z^ 



(A13) 
(A14) 
(A15) 



as well as the same terms with {ijk) permuted. However, 
it is easy to see that the second two forms reduce to the 
the first since expression IA14I is equivalent to 



-1 



-1 



(zi - z J ) 2 {z j - z k ) (zi - Zj)(zj - z k ) 2 
and expression I A 1 51 is equivalent to 
-1 -1 



(Z k - Z l ) 2 (z j - Z k ) (Zi ~ Zj)( Zk - Zi) 2 



(A16) 



(A17) 



It follows then, that it's enough to consider forms of type 
of Eq. IA13I for triplets (ijk) (as well as the same form 
with permutations of (ijk)). When (ijk) — > z, the most 
general zero energy eigenstate should have the form 



<t>G <x 



(Zi - Zj) 2 ( Z j - Z k ) 



(A18) 



(or a form like this with any permutation of (ijk)) where 
F(. . .) must be analytic (i.e., with no poles). 

Now arbitrarily pair up and relabel the particles, e.g. 
(za 1 ,z Bi ), (za 2 ,zb 2 ), (z An/2 ,z Bn/2 )- Look at the 
most singular part of </>g as particles within these pairs 
approach each other, while pairs are kept separated. 



1 



1 



t>a oc 



(z Al - Z Bl ) 2 (ZA 2 - Z B2 ) 2 

1 



(z Ak _ - z Bk _Y 



(A19) 



Since we've isolated the most singular part of </>g, it's 
clear that </>w cannot contain factors (z^ — z Bi )~ x ■ If 
we now consider triplets (Ai,Bi,k) Vk ^ {Ai,Bi}, and 
bring particle k close to (A\, B\) it's clear that <f>N_ must 

contain a factor of (za x — z k )~ x or (z Bl — Zk)^ 1 , but not 
both, in order to satisfy the requirements on the pole 
structure deduced above. We might be led to naively 
define 



Em 
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(zA t - z v j )(z c j - z Bl ) 



(A20) 



where C j U V' J = {zj; C j H V j = 0, (i.e. a partition of 
the set of particle coordinates), and j indexes all possible 
partitions. However, the pole structure places further 
restrictions on the sets C and T>. In particular, Ai and 
Bi cannot both be in C or in V. Otherwise, supposing 
Ai,Bi £ C, we have (schematically) the following: 



(f> G OC 



(z Al - z Bl ) 2 " ' (z At - z Bi ) 2 

1 

' (ZAi ~ Z Ai )(z Al - ZBi) " 



(A21) 



and we immediately recognize, that the triplet 
(Ai,Ai, Bi) has too many poles (Q < —3). We conclude, 
that for i-th pair only one factor is allowed in (f> « , either 

(za 1 — Z Ai ) , or (z Al — z Bi )~ x . That is, the partitions 
are such that C includes only one representative of any 
pair, and T> includes the complementary member of this 
pair, i.e. necessarily A; e C, Bi G V or Bi e C, Ai e V. 
At this point we can recognize that our notation of C and 
T> is redundant, and that we can rewrite 



>n oc 



E 

fc£Partitions 



n 

&3 



4>k 



(z A * 



(A22) 



where <p k cannot contain any more poles, and the conven- 
tions are that for all k the same coordinates are paired 
up, i.e. {A^,B^} = {Af,Bf'} are equal as sets. The 
difference between partition k and partition k' is the 
order of coordinates within a pair, i.e. we could have 
z A k — Zgk' , z B k = z j^k' , or Zj^k — z j^y , z B k — Zgk' • 
Clearly the sum over k is a subset of symmetrization 
over all particles. 

Finally, we should also include the exchange of pairs 
(z A k, z B k) <-> (z A k, z B k) 7 since the most singular part is 

i i j j 

symmetric under this exchange, and arrive at 
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e En 



pairings /cGPartitions 

Then for the particular choice of pairs we'll have 



4>k 

(z A k — Z B k) 



(A23) 



e e n 

pairings fcGPartitions 



0k 



(z A k — Z B k) 



(A24) 



and recover the whole eigenfunction by symmetrization over all particles and multiplication by the Jastrow factor 
squared: 



e e n 



(z A k — Z B k) 



(A25) 



pairings fcgPartitions 

The wavefunction of the densest state has as few zeros as possible, and to find it we may choose (j>k = 1, then 



I pairings fcePartitions \ i K A i B i ' i=£j V A i B j ' J 



n i n - 



(A26) 
(A27) 



This is can be recognized as the proposed Gaffnian wave- 
function with no broken pairs and no added flux. To 
obtain the states of lower density we need to consider 
the case of non-constant 4> k (z Al , z Bl ; . . . ; z An/2 , z Bn/2 ). 
By analyzing the symmetry of denominators we find 
that <pk must be symmetric under the exchange of pairs 
(z Ai ,z Bi ) «-> (z Aj ,z Bj ). We now claim that a complete 
basis for functions that satisfy this symmetry condition 
is given by functions of the form 

N/2 

E tlM"A Tm ,z BrW ) (A28) 

rGSjv i=l 

2 

where the fiS are chosen from a basis for arbitrary poly- 
nomials of their two arguments. While this may seem 
to be a strange way to write a basis for the polynomial 
4>k (z Al , zbx ; ■ ■ ■ ; za n/2 , zb n/2 ) this is actually a form well 
known to physicists. To sec this, imagine a system of N/2 
bosons where the "position" of each boson is specified by 
two coordinates (z±, z 2 ). The functions are basis func- 
tions for the single "particle" positions. All multiparticle 
states can be written as symmetrized (bosonic) linear 
combinations of the occupancies of these basis states. 

Consider now the case, when we've added n quanta of 
flux to the Gaffnian ground state. The highest degree of 
fi() is n, and we could choose basis polynomials fi(z\, z 2 ) 
of the form z" 1 ^ 12 with < ni,n 2 < n. The dimension 
of this space is (n + l) 2 . However, a different basis set 



turns out to be more useful. Specifically, it is useful to 
separate functions ft that vanish in the limit Z\ — » z%, 
from ones that do not. 

We choose a basis for our space of /j which decom- 
poses into two disjoint sub-bases: the symmetric z™ 1 2 + 
z" 2 ^ 1 with < rii < ri2 < ti and the antisymmetric 
z" x z™ 2 - z^z^ 1 with 0<7ii<n2<n-l. The dimen- 
sions of subspaces spanned by them are h(n + 2)(n + 1) 
and \ {n + l)n respectively. Clearly the span of the anti- 
symmetric sub-basis vanishes as z\ — > z 2 . The quotient of 
the full space by the span of the antisymmetric sub-basis 
is just the span of the symmetric sub-basis (S), i.e. sym- 
metric polynomials. Of these, polynomials which vanish 
as z\ — ► Z2 are spanned by [z\ — z 2 ) 2 (z™ 1 z 2 2 + z™ 5 ^ 1 ), 
with < rii < ri2 < n — 2. The dimension of this sub- 
space is hn(n — 1). 

The quotient of <S by the subspace of the vanishing 
symmetric polynomials has dimension 2n + 1 and con- 
tains symmetric polynomials in 2 variables that do not 
vanish in the limit z\ — > z 2 , we'll call this quotient Q. 
However, by considering the Taylor expansion^* of Read- 
Rezayi pairing form $(z 1; z 2 ) given in Eq. IA2I we've al- 
ready found a set of 2n + 1 linearly independent symmet- 
ric polynomials in 2 variables, thus we may choose them 
as the basis of this quotient space. 

Now given a choice of 4>k we obtain a zero energy state 
of the Hamiltonian. Further, all possible zero energy 
states can be written in this way. We can now decompose 
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any into basis polynomials fi of the above described 
form. Let our choice be such that /*() for 1 < i < -j 
belong to the subspace of polynomials that vanish as 
z\ — » zi and for €■ + 1 < i < y belong to the com- 
plementary subspace, i.e. Q. Then each vanishing /;() 
simplifies with the appropriate factor in the denomina- 
tor of 4>g producing a "broken pair" , and the remaining 
factors form what we above called 0(), whereas the prod- 
uct of non- vanishing /j() can be reexpressed as a linear 
combination of Read-Rezayi pairing forms <!>(). Thus we 
conclude that the most general zero energy eigenstate of 
Hq is of the conjectured form, and therefore we counted 
the complete degeneracy of eigenstates for a given value 
of additional flux n. 



APPENDIX B: AN EXAMPLE OF ANGULAR 
MOMENTUM ADDITION 

We would like to determine the full angular momentum 
spectrum of the zero energy states of the Hamiltonian 
Pg using Eq. [SJ Here we will consider the example of 
N = 4 particles and n = 3 (6 quasiholes). Eq. [5] tells 
us that we should have a total number of zero energy 
states given by the sum of three terms corresponding to 
F = 0, 2,4. For F — 4 we have (6 bosons in 1 orbitals) 
® (4 fermions in 4 orbitals). Both 6 bosons in 1 orbital 
on 4 fermions in 4 orbitals have L = 0, so overall this is 
an L = state. The F — 2 case is more tricky. Here we 
have (6 bosons in 2 orbitals) ® (2 fermions in 3 orbitals). 
First we take 6 bosons in 2 orbitals. When there are two 
orbitals on a sphere, we clearly have L = 1/2. So the 
two orbitals have L z = ±1/2. There are 7 ways to fill 



these orbitals with 6 bosons, which we can write as (6,0); 
(5,1), . . . (0,6). Counting the total L z of each of these 
states, we get 3, 2, 1, 0, —1, —2, —3 which we recognize as 
being L — 3. Thus, 6 bosons in 2 orbitals is L = 3. 
Similarly, we look at 2 fermions in 3 orbitals. The three 
orbitals must be L = 1 with L z — 1, 0, — 1. We can fill the 
three orbitals with two fermions in 3 ways, which have 
L z = 1,0,-1 so we recognize this as L = 1. Now we 
must add together the angular momentum of (6 bosons 
in 2 orbitals) ® (2 fermions in 3 orbitals). This means 
we need to add L = 3 with L = 1. By the usual angular 
momentum addition rules we obtain L — 2,3, 4. Finally, 
we turn to the F = case. Here we have (6 bosons in 
3 orbitals) <g> (0 fermions in 2 orbitals). The fermions 
in 2 orbitals clearly has L — 0. It is a simple exercise 
to count up the possibilities for 6 bosons in 3 orbitals. 
We discover that this has L — 0, 2, 4, 6. Putting together 
all of the results we find that the zero energy states of 
the Hamiltonian Pg for N particles with n — 3 occur at 
angular momentum L = 0,0,2,2,3,4,4,6 which agrees 
with the results of exact diagonalizations. 

APPENDIX C: FURTHER GENERALIZED 
WAVEFUNCTIONS 



Although there may be many possible ways to gener- 
alize Gaffnian wavefunction 18 , the form written in Eq. ^ 
suggests a generalization from paired to clustered wave- 
functions where instead of dividing the particles into two 
groups, we divide the particles into g-groups. Let us as- 
sume the number of particles N in the system is divisible 
by g and write N = gn. We then write the wavefunction 



* = S 



n 

a=l 



n 



(a— l)n<i< < 7 <an 



n 



l<a<b<g 



n 



_^ z (a— l)n+i z [b~ l)n-\-i 



n 



(Cl) 



l<i<j<N 



with m > 1 where again S symmetrizes or antisym- 
metrizes for bosons (odd m) or fermions (even m) re- 
spectively. Counting powers of z we discover that this 
wavefunction occurs at flux 

iV = (N/g-l)-(g-l)+m(N -1) (C2) 
= (l/g + m)N-(g + m) (C3) 

corresponding to a filling fraction v = g/(gm + 1) which 
is just the Jain sequence. Furthermore the precise value 
of the flux (the shift) is also in agreement with the Jain 
series. This construction clearly reproduces the Gaffnian 
for g = 2. For the bosonic case (m = 1) for general 
g this construction produces a wavefunction that does 



not vanish when g particles come to the same point, but 
vanishes as g + 1 powers as the g + l'st particle arrives 
at that point. However, for g > 2 this trial wavefunc- 
tion is not the densest possible wavefunction with this 
particular property^. Nonetheless, we believe that this, 
and other related wavefunctions can generally be con- 
structed with simple projection rules. For example, for 
the g = 3, m — 1 case of Eq. ICll this wavefunction is the 
unique densest wavefunction that does not vanish as 3 
particles come together, that always vanishes as at least 
4 powers when 4 particles come together, and vanishes 
faster than 4 powers if particles are brought together in 
groups of 2 and then two groups of 2 are brought to- 
gether. 
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